/iD-ni23  986  STRESS  HERSUREHENTS  IN  IC£(U>  COLD  REGIONS  RESEARCH  AND 
ENGINEERING  LAB  HANOVER  NH  G  F  COX  ET  AL.  AUG  83 
CRREL-83-23 


1/1 


UNCLASSIFIED 


F/G  8/12  NL 


For  convers/on  of  SI  metric  units  to  U.SJBritish 
customary  units  of  measurements  consult  ASTM 
Standarri  B380,  Metric  Practice  Quide,  published 
by  the  American  Society  tor  Testing  and  Materi¬ 
als,  1916  Bace  St,  Philadelphia,  Pa.  19103. 


Cover:  Cylindrical  steel  Ice  stress  sensor. 


CRREL  Report  83-23 

August  1983 


Stress  measurements  in  ice 


Gordon  F.N.  Cox  and  Jerome  B.  Johnson 


Accession  For 


NTIS  GRA&I 
DTIC  TAB 
Unannounced 
Justif  ication_ 


n 


By. - 

Distribution/ 

Availability  Codes 
|Avail  and/or 
Special 


lor  publle  nrioMO;  dlotrtbutlon  unitmitod 


SBCUHITV  CLAmnCATION  or  THIS  PAGE  fWhan  OaM  Bntta4> 


REPORT  DOCUMENTATION  PAGE 


READ  INSTRUCTIONS 
BEFORE  COMPLETING  FORM 


AiTFITlf 


CRKEL  Report  83-23 


4.  TITLE  raiE  SvMa<t) 

STRESS  MEASURHIIENTS  IN  ICE 


AUTHONfA) 


Gordon  F  Jf.  Cox  and  Jerome  B,  Johnson 


psiiromuNa  organization  name  and  address 

UiS.  Army  Cold  Regions  Research  and  Engineering  Laboratory 
Hanover,  New  Hampdiire  037SS 


It.  controlling  office  name  and  address 


S.  TYPE  OF  REPORT  A  PERIOD  COVERED 


s.  performing  org.  report  NUMGER 


8.  CONTRACT  OR  GRANT  NUMBERfa) 

U.S.D.I  Order  No.  2LA6000-1088 


to.  PROGRAM  ELEMENT,  PROJECT,  TASK 
AREA  A  WORK  UNIT  NUMBERS 

lUR  4A161 10li491D.  Task  00, 

Work  Unit  388 


12.  REPORT  DATE 


U,5.  Amqf  Odd  Regions  Researdi  and  Minerals  Management  Service  August  1983 _ 

EngpwAAring  Inboratoiy  and  U.S.  Department  of  Interior  o-  number  of  pages 

Hmover,  New  Hampdiire  037SS _ Reston,  AArginia  22091 _ _ 


CV  NAME  A  ADORESarr/  dlUtimt  inm  CanNofflna  ONleaJ  IS.  SECURITY  CLASS,  (ol  M»  twport) 

Unclassified 


mm 


sTaI  i7Tn»rj*l 


nnUTlON  STATEMENT  (ml  Mm  Hmfort) 


Appnmd  for  puUic  release;  distribution  unlimited. 


tr.  MSTRIEUTION  STATEMENT  (ml  IM  < 


I  aaMaaW  la  WoeA  SO,  II  dKfatanl  Aeoi  Rmpott) 


W.  KEYWONDGrCi 
Cold  regions 
be 

loe  stress 
Laboratory  tests 
Stress  sensors 


fCawfIniia  on  raaoroa  mUm  linmmmmmarr  ooF  IBMNiy  Ap  MoeA  nuaSorJ 


1  i  1 »  7,1  “XM  Z 


”~The  proUems  associated  with  measuring  stresses  in  ice  are  reviewed.  Theory  and  laboratory  test  results  are  then  pre¬ 
sented  for  a  stiff  c)dindrical  semor  made  of  steel  that  is  desipied  to  measure  ice  stresses  in  a  biaxial  stress  field.  Load¬ 
ing  testa  on  fteahwater  and  saline  ice  blocks  containing  the  biaxial  ice  stress  sensor  indicate  that  the  sensor  has  a  resolu¬ 
tion  of  20  kft  and  an  accuracy  of  better  than  15%  under  a  variety  of  uniaxial  and  biaxial  loading  conditions.  Principal 
straaa  dbections  can  also  be  determined  witiiin  y.  The  biaxial  ice  stress  sensor  is  not  significantly  affected  by  variations 
In  the  ice  elastic  modulus,  ice  creep  or  differential  thermal  expansion  between  the  ice  and  gauge.  The  sensor  also  has  a 
low  temperature  sensitivity  (S  kPi/^. 

-  ? 


BB.STnMn 


<SF  t  Noveeiai 


Unclassified _ 

SECURITY  CLASSIFICATION  OF  THIS  PAGE 


flRMn  Data  Bnlmrmd) 


PREFACE 


This  report  was  prepared  by  Dr.  Gordon  F.N.  Cox,  Geophysicist,  of  the  Snow  and  Ice  Branch, 
Research  Division,  U.S.  Army  Cold  Regions  Research  and  Engineering  Laboratory  and  Dr.  J.B. 
Johnson,  Geofdiysicist,  Geophysical  Institute,  University  of  Alaska.  The  work  was  jointly  sup¬ 
ported  by  the  Minerals  Management  Service  of  the  U.S.  Department  of  the  Interior,  Order  No. 
2LA6000-1 088,  Development  and  Testing  of  a  Field  Ice  Stress  Measurement  System,  and  the 
UJS,  Army  Cold  Regions  Research  and  Engineering  Laboratory,  In-House  Laboratory  Indepen¬ 
dent  Research,  ILIR,  DA  Project  4A1 61 101 A91 D,  Task  00,  Work  Unit  388,  Ice  Stress  Meter. 

The  authors  appreciate  the  assistance  provided  by  Nancy  Perron,  Steve  Decato  and  Bill  Bos- 
worth  in  growing  and  handling  the  large  ice  blocks,  as  well  as  preparing  ice  thin  sections.  Larry 
Gould  designed  the  biaxial  loading  machine  which  was  then  fabricated  by  Bill  Burch.  The  authors 
also  thank  Dr.  Dev  Sodhi  of  CRREL  and  Lh.  Ryush  Dutta  of  IRAD  Gage  for  technically  reviewing 
the  manuscript  of  this  report. 

The  Contents  of  this  repwt  are  not  to  be  used  for  advertising  or  promotional  purposes.  Citation 
of  brand  names  does  not  constitute  an  official  endorsement  or  approval  of  the  use  of  such  com¬ 
mercial  products. 


CfWnCNTS 


Page 

Abstract .  i 

Rtefoce .  ii 

Introduction .  1 

ftevious  work .  1 

Stress  measunments .  1 

Design  considerations .  2 

Stress  sensors .  2 

Biaxial  ice  stress  sensor . 5 

Biaxial  stress  sensor  theory .  6 

Gauge  deformatkHi .  6 

Stresses  associated  with  cylindrical  sensors .  7 

Determirution  of  ice  stresses .  1 2 

Gauge  calibration .  14 

*lvahiatioo  of  the  biaxial  ice  stress  sensor . . .  1 S 

Temperature  sensitivity .  16 

Biaxial  loading  test  equipment .  16 

Biaxial  loading  test  results .  21 

DifTerential  thermal  expansion .  28 

Longterm  drift .  28 

Diacuiaion  of  test  results .  29 

Omdusions .  29 

literature  cited .  30 

ILLUSniATIONS 

Rgure 

1.  Haxial  ice  stress  sensor .  5 

2.  Scfaenutic  of  biaxial  ice  stress  sensor .  6 

3.  Ran  view  of  cylindrical  sensor  frozen  into  ice. .  8 

4.  Nomnlized  radial  stress  at  ice-gauge  boundary  parallel  to  uniaxial  loading  direction 

versus  ice-gauge  modulus  ratio .  10 

5.  Normalized  stress  distribution  in  the  ice  surrounding  the  biaxial  ice  stress  sensor 

under  uniaxial  loads .  1 1 

6.  Nomulized  radial,  tangential  and  diear  stress  distribution  in  the  ice  surrounding  the 

biaxial  ice  stress  sensor  under  biaxial  loads .  12 

7.  Normalized  radial,  tangential  and  drear  stress  distribution  in  the  ice  at  the  ice-gauge 

boundary  under  uniaxial  and  biaxial  loads .  12 

8.  Variation  of  A  with  ice  moduli .  14 

9.  Variation  of  A  with  ice  moduli .  14 

10.  ffydrauBc  pressure  cell  used  to  calibrate  the  biaxial  ice  stress  sensor .  IS 

1 1 .  Variation  of  wire  pniod  with  temperature  for  each  of  the  three  wires  in  the  biaxial 

ke  stress  sensor .  16 

12.  Freezing  chamber  used  to  grow  ice  blocks  for  the  stress  sensor  verification  tests .  1 7 

13.  Qose>up  of  fieezing  diamber  showing  aluminum  coldplate  and  ice  block .  17 

iii 


Figure  Page 

14.  Vertical  thin  section  from  fierh  water  ice  block .  18 

15.  Support  frame  used  to  position  the  sensor  in  the  hole  while  it  was  frozen  in  the  ice 

block .  19 

16.  ffiaxial  loading  machine  used  in  stress  sensor  verification  tests .  19 

17.  Ice  block  and  sensor  in  biaxial  loading  machine .  20 

18.  Equipment  used  to  measure  block  strains  during  testing .  20 

19.  Position  of  sensor  relative  to  the  loading  directions  for  each  of  the  four  ice  blocks 

tested .  21 

20.  Measured  stress  versus  applied  stress  and  block  strain  for  block  1 .  22 

21 .  Measured  versus  applied  stress  in  the  i4  and  B  directions  for  block  2 .  23 

22.  Measured  versus  applied  stress  for  block  3 .  25 

23.  Measured  versus  applied  stress  for  block  4 .  28 

24.  Variation  of  wire  period  with  time  for  each  of  the  three  wires  in  the  biaxial  ice 

stress  sensor .  28 


STRESS  MEASUREMENTS  IN  ICE 


Gordon  F.N.  Cox  and  Jerome  B.  Johnson 


INTRODUCnm 

Relud>le,  inexpensive  ice  stress  measurements  are  needed  to  solve  a  variety  of  ice  related  prob¬ 
lems.  These  include:  measuring  and  monitoring  ice  loads  on  marine  and  hydraulic  structures; 
determining  the  magnitude  of  ice  forces  associated  with  ice  drift,  ride-up,  pile-up  and  pressure 
ridge  formation;  measuring  thermal  ice  pressures  in  reservoirs;  and  assessing  the  effects  of  ice  con¬ 
vergence  on  the  performance  of  large  icebreakers  and  tankers. 

Researchers  have  obtained  estimates  of  ice  loads  on  structures  by  considering  the  failure  strength 
of  the  ice;  however,  because  of  the  uncertainty  in  the  large  scale  mechanical  properties  of  the  ice 
sheet,  these  estimates  may  be  too  conservative.  In  situ  measurements  of  stress  in  ice  are  needed  to 
accurately  determine  ice  loads  on  structures. 

In  this  report  we  first  review  the  problems  of  measuring  ice  stress,  together  with  the  findings 
and  accomfdishments  of  other  investigators  who  have  worked  on  the  development  of  ice  stress 
senson.  We  then  present  theory  and  laboratory  test  results  for  a  stiff  steel  cylindrical  sensor  de- 
s^ned  to  measure  ice  stresses  in  a  biaxial  stress  field. 


PREVIOUS  WMK 
StKSi  measnrancnti 

Hie  general  problems  of  measuring  stresses  in  ice  and  other  materials  have  been  adequately 
addressed  by  Metge  et  al.  (1975).  Stress  in  any  material  cannot  be  measured  directly.  It  must  be 
determined  by  measuring  the  strain  deformation  of  the  material  or  by  measuring  the  strain  of  an 
elastic  inclusion  embedded  in  the  medium.  In  materials  that  deform  elastically  and  where  the  elastic 
modulus  is  known,  the  stress  can  be  calculated  given  the  material  deformation  and  elastic  modulus. 
However,  in  materials  such  as  ice,  which  exhibit  time-dependent  deformation  and  wide  variations 
in  the  elastic  modulus,  an  imbedded  inclusimi  must  be  used  to  measure  the  stress  where  the  stress- 
strain  relationship  and  the  inclusion  factor  of  the  sensor  are  known.  The  inclusion  factor  is  de¬ 
fined  u  the  ratio  of  the  undisturbed  ice  pressure  to  the  pressure  felt  by  the  sensor.  It  should  be 
noted  that  the  elastic  modulus  of  ice  can  vary  over  one  order  of  magnitude,  from  0.7  to  10  GPa, 
depending  on  the  ice  salinity,  temperature,  grain  size,  crystal  orientation  and  strain  or  loading 
rate  (Traetteberget  al.  1975,  Vaudrey  1977,  and  Schwarz  and  Weeks  1977). 


If  the  ebstic  modulus  of  the  inclusion  or  stress  sensor  is  different  than  that  of  the  surrounding 
ice,  the  sensor  will  change  the  local  stress  field  in  the  ice.  A  sensor  that  is  stiffer  than  the  ice  will 
support  some  of  the  load  that  would  otherwise  be  supported  by  the  surrounding  ice.  It  will  con¬ 
centrate  the  stress.  If  it  is  softer  than  the  ice  it  will  deflect  easily,  requiring  the  sunounding  ice  to 
support  more  of  the  load.  In  addition  to  the  rebtive  stiffness  between  the  sensor  and  the  ice,  the 
inclusion  factor  of  the  sensor  depends  on  the  inclusion  geometry  and  direction  of  the  applied  stress 
in  the  ice.  Hence,  precise  knowledge  of  the  strain  in  the  sensor,  the  stress-strain  relationship  of  the 
sensor  material  and  the  inclusion  factor  are  needed  to  obtain  the  stress  in  the  ice.  The  main  problem 
b  to  determine  the  inclusion  factor  and  to  design  the  sensor  in  such  a  way  that  it  remains  nearly 
constant,  even  if  the  ice  properties  change.  In  the  design  of  ice  stress  sensors,  variations  in  ice 
modulus  can  be  dealt  with  by  either  selection  of  a  thin,  wide  sensor  (Metge  et  al.  1975,  Templeton 
1979)  or  a  sensor  that  is  much  stiffer  than  the  ice  (Nelson  et  al.  1977,  Johnson  and  Cox  1980). 

Other  potential  probbms  in  measuring  ice  stress,  dted  in  part  by  Metge  et  al.  (1975)  and  ebb- 
orated  by  Templeton  (1979),  include  the  effects  of  nonebstic  behavior,  differentbl  thermal  ex¬ 
pansion  between  the  ice  and  the  gauge,  and  overloading  of  the  surrounding  ice.  Two  examples  of 
nonebstic  behavior  include  ice  creep  and  localized  plastic  yielding  around  the  sensor.  Creep  tends 
to  reduce  the  ebstic  modulus  of  the  ice  and  to  increase  its  varbbility.  Localized  plastic  yielding 
around  the  sensor  may  produce  deformation  in  the  sensor  that  remains  even  after  the  applied  or 
far-fkld  stress  is  removed.  Differential  thermal  expansion  between  the  ice  and  sensor  takes  pbce 
whenever  the  system  temperature  changes.  This  may  result  in  anomalous  stress  readings  entirely 
caused  by  the  different  thermal  expansion  characteristics  between  the  ice  and  gauge.  If  the  sensor 
produces  a  hi^  stress  concentration  in  the  surrounding  ice,  overloading  of  the  ice  may  result  in 
localized  failure  of  the  ice  around  the  sensor  and  further  stress  measurement  errors.  Templeton 
(1979)  suggesb  that  all  these  probbms  can  be  minimized  by  the  choice  of  a  thin,  wide  sensor  hav¬ 
ing  a  modulus  close  to  that  of  ice.  Johnson  and  Cox  (1980)  also  demonstrated  that  these  problems 
are  not  apparent  for  stiH  cylindrical  sensors. 

Oeatpi  oonsUeiatioiia 

b  l^t  of  die  previous  discussion  and  recommendations  by  Templeton  (1979),  the  following 
are  important  design  considerations  for  an  ice  stress  sensor: 

1 .  Hie  sensor  should  not  be  affected  by  varbtions  in  the  ice  ebstic  modulus  and  by  nonebstic 
behavior  of  the  ice. 

2.  Hie  sensor  should  have  a  low  temperature  sensitivity  and  not  be  significantly  affected  by 
differential  thermal  expansion. 

3.  Hie  wnsor  diould  not  gready  overload  the  ice. 

4.  The  sensor  diould  be  appropriately  sized  and  be  rugged  and  leakproof. 

5.  Hie  sensor  diould  be  inexpensive,  easily  installed  and  monitored,  and  have  a  stable,  repeat- 
abb  response. 

Numerous  attempb  have  been  made  to  derign  such  a  sensor.  They  are  described  below. 


Sbcas  season 

ke  stress  sensors  have  been  developed  by  Esso  Resources  Canada  (formerly  Imperial  Oil  Limited, 
K)L),  the  University  of  Alaska,  Exxon  IVoduction  Research  (EPR),  the  National  Research  Council 
of  Canada  (NRC)  and  Oceanographic  Services  Inc.  (OSO  in  cooperation  with  IRAD  Gage.  Hawkes 
(1969b)  also  used  a  idiotoebstic  streasmeter  to  measure  stresses  in  frozen  sands  and  Baumarm 
(1979)  used  earth  pressure  celb  to  measure  stresses  in  river  ice.  These  sensors  vary  widely  in  geom¬ 
etry  and  modulus. 

The  lOL  sensor  desoibed  by  Metge  et  al.  (197S)  b  a  thin,  wide,  soft  sensor  having  an  effective 
ebstic  modulus  bas  duui  that  of  ice.  It  conrists  of  a  double  sandwich  of  aluminum  pbtes  and 
elastomeric  material  dut  deforms  under  ai^lied  stress.  The  amount  of  deformation  is  determined 
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^  OMlweBMOt  of  the  change  in  capacitance  between  the  metal  plates.  The  gauge  is  0.79  cm 
Aide.  122  cm  wide  and  is  designed  to  extend  duough  the  full  thickness  of  an  ice  sheet.  The 
pngs  has  been  wide^  used  by  lOL  to  meanire  the  ice  stress  around  man^nade  fill  islands  in  Mac- 

KemdeBdy* 

Ihe  IMmiaity  of  Aladca  puge  is  described  by  Nelscm  et  al.  (1977).  It  consists  of  a  2.S4-cm- 
diamettr  dundnun  cylhider  that  is  7.62  cm  long,  with  a  S.08-cm4ong  by  1 .27-cm-diameter  re¬ 
duced  aeelioo.  Fbur  strain  puges,  parallel  and  perpendicular  to  the  axis  of  the  cylinder,  are  con¬ 
nected  ha  a  hrUp  to  read  toisioii  or  compreadon  in  the  bar  and  provide  temperature  compensa- 
tioB.  The  ahmdnura  cylinder  fits  inside  a  copper  tube  and  has  1 27-cm-diameter  steel  bolts  on 
each  end  to  grip  the  ice.  The  eoppet  tube  is  sealed  to  the  aluminum  bar  with  silicone  rubber  and 
the  entire  asaonhly  is  coated  with  silicone  rubber.  Both  aluminum  and  brass  have  been  used  to 
eoBStmet  the  p«^.  The  dimensions  of  the  pup  have  also  been  sh^tfiy  modified  in  various  field 
measuteaiMit  pogtams. 

The  Unhferdty  of  Alaska  pup  is  frozen  into  the  ice  horizontally  and  stresses  are  measured  at 
the  ends  of  the  qrhnder.  It  is  a  uniaxial  device,  in  that  it  can  only  accurately  measure  ice  stress 
when  ortented  parallel  to  a  uniaxial  stress  field.  In  sudi  an  orientation  the  pup  has  an  averap 
inchiaioB  helm  of  1 :32.  Other  tests  have  shown  the  sensor  to  have  a  transverse  sensitivity  of  up 
to  25%  of  the  applied  load.  We  should  point  out  that  because  the  off-axis  loading  characteristics 
ef  file  gaup  have  not  yet  been  fully  evaluated,  measuremoits  in  biaxial  stress  fields  only  provide 
an  hidintion  of  the  state  of  stress,  not  values  of  fiie  magnitude  and  direction  of  the  principal 
stremes.  in  uniaxial  stress  fields,  when  the  pup  is  oriented  at  some  angle  to  the  applied  stress, 
fip  nsultsalao  have  to  be  carefully  interpreted.  Nevertheless,  the  pup  has  been  successfully 
used  to  understand  the  genoal  magnitude  of  the  stresses  in  ice  around  man-made  and  natural 

structures  (SackiBgR  and  Nelson  1979a,  b.  Nelson  and  Saddnpr  1976). 

Use  EPR  ice  pressure  sensor  consists  of  a  thhi,  wide  panel  fiiat  extends  through  the  full  fiiick- 
ness  of  tte  ice  fiwet  (Temfdeton  1979).  It  is  1.11  cm  fiiick  and  about  45  cm  wide.  Aluminum 
season  wifii  strahi  pages  ate  used  to  measure  ice  pressures  normal  to  the  panel.  Constructed 
season  have  exhMted  en  effective  elastic  modulus  of  (9  to  1 .84  GBa,  close  to  that  of  sea  ice. 

Numerous  HUi^ytical  studies  have  been  prformed  to  determine  fiie  indusion  factor,  transverse 
sanalfiviw  and  diffnentid  thermal  expandon  characteristics  of  the  EPR  pressure  sensm.  Chen 
(1981a)  fi»vn  that  the  aveiapinduakm  factor  of  the  pup  is  dose  to  1:1.  Chen’s  elastic  finite 
slaaiaBtaaa^yasaabo  fiiow  fiiat  the  pi^  produces  a  nuudmum  stress  concentrafim  factor  of  1.S 
in  file  ioa  near  fiie  panel’s  edp.  Other  ftilte  elemoit  amdyses  performed  by  Qien  (1981b),  using 
aa  dnilc  gtetic  ke  model,  indicate  fiiat  the  sensor  hu  a  very  low  transverse  sensitivity.  Trans- 
vane  prsaiaiaa  ip  to  1^  MPs  result  hi  maximum  anoimlous  pressure  on  the  sensing  face  of  only 
0.17  Iffh.  He  found  transverse  presiure  effects  to  be  greatest  at  small  transverse  pressures. 

hi  addifiOB  to  fiieae  studies,  Templeton  (1981)  conducted  an  elasto^tatic  analysis  to  determine 
file  gpup’s  seaaitMty  to  differentid  thermal  erqianskm  between  the  ice  and  pup.  His  results  show 
thatenoR  fiom  diffeieafial  fiiermal  expansion  are  independent  of  the  senior’s  elastic  modulus  and 
can  be  mjahahad  by  dioosii^  a  fiiin,  wide  sensor.  He  condudes  fiut  errors  due  to  transverse 
praanires  and  differential  thermal  expansion  total  leas  than  10%  of  the  measurement. 

Chen  (1981a)  and  Chen  and  Tmipleton  (1983)  also  present  some  results  of  field  verification 
tests  OB  file  EPRioe  preanue  tensor.  In  these  tests  fiie  sensor  was  frozen  into  larp  sea  ice  blocks 
meawaiag  3  by  6  m  by  the  fiill  fiiicknesi  of  the  ice  fiieet.  The  Mocks  were  loaded  using  two0.4S- 
HN  caparftylqfdisidic  qdiiiders.  The  test  results  showed  that  the  measured  stress  was  within  15% 
of  the  Bonfiaal  antic^arted  response  for  aU  apjdied  ice  pressures  greater  than  0.69  MPa.  At  pres- 
sans  lawfiimi  0j69  MFh,  the  actual  response  ivf  the  pup  was  much  less  than  predicted.  No  field 
teetehan  bean  made  hi  vnif  he  analytic*’  adin  of  the  sensor’s  transverse  sensitivity. 

Barth  pnmae  edbhava  dk.  '  at  o  moaftm  stresses  in  ice.  Bautnaim  (1979)  used  Terra- 
Tiahnoinp  earth  pressure  cells  w  r.ieesure  ice  stremes  in  the  St.  Marys  River,  Midiigan.  Theoells 
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coimt  of  two  sealed  20-  by  30-  by  l-cm  steel  plates  filled  with  hydraulic  fluid  which,  in  turn,  is 
pressurized  by  CO]  gas.  Pressures  are  monitored  with  a  gas  regulator.  The  sensors  have  a  sensi¬ 
tivity  of  690  Ph  and  an  operating  range  of  0  to  1 .36  MPa.  Baumann  calibrated  the  cells  in  the 
laboratory;  however,  he  does  not  give  details  of  the  calibration  procedure.  As  no  attempt  was 
made  to  determine  the  transverse  sensitivity  of  the  cells  in  ice,  measurements  of  stress  need  to  be 
carefully  interpreted,  as  with  the  University  of  Alaska  gauge. 

The  NRC  gauge  is  a  thin-walled  aluminum  tube  having  an  outside  diameter  of  S.0  cm,  a  wall 
thickness  of  0.3  cm  and  a  length  of  10.0  cm.  Three  strain  gauges  are  circumferentially  bonded  to 
the  inside  of  the  tube  at  120°  intervals  to  measure  the  tube  deformation  under  stress.  The  stresses 
on  the  tube  are  calculated  from  the  circumferential  strains,  assuming  that  the  ice  behaves  as  a  lin¬ 
ear  elastic  material  (Frederking  1980). 

The  tube  material  and  dimensions  are  chosen  to  minimize  the  stress  and  strain  distributions  in 
the  surrounding  ice  associated  with  the  presence  of  the  sensor.  The  displacements  of  the  outer 
diameter  of  the  tube  approximate  those  of  a  solid  ice  cylinder  of  the  same  size.  In  effect,  the 
gauge  has  an  inclusion  factor  close  to  1 : 1 . 

Both  laboratory  and  field  tests  of  the  NRC  puge  show  that  a  tubular  transducer  can  success¬ 
fully  be  used  to  determine  the  principal  stress  direction.  This  agrees  with  the  theory  presented 
by  Frederking.  fat  calibration  tests  where  the  elastic  modulus  of  the  ice  was  known,  the  measured 
principal  stresses  were  within  20%  of  the  applied  stress.  In  the  field,  poorer  agreement  was  ob¬ 
tained;  however,  this  was  in  part  due  to  the  uncertainty  in  the  applied  stress  field. 

Calculations  presented  later  in  this  report  show  that  the  inclusion  factor  of  a  cylindrical  inclu¬ 
sion,  having  an  effective  modulus  close  to  that  of  ice,  is  very  sensitive  to  small  changes  in  the  ice 
modulus.  If  the  ice  creeps,  resulting  in  a  large  decrease  in  the  effective  ice  modulus,  significant 
enors  can  be  expected  with  this  type  of  sensor. 

The  OSI  sensor  described  by  Johnson  and  Cox  (1980)  is  similar  to  the  NRC  puge  in  that  it 
also  is  a  cylindrical  inclusion.  The  sensor  has  an  outer  diameter  of  2.86  cm,  a  wall  thickness  of 
0.79  cm  and  a  length  of  57.0  cm.  The  ends  of  the  pup  are  fitted  with  rounded  end  caps.  Rela¬ 
tive  to  the  NRC  puge,  this  sensor  is  much  stiffer,  having  an  effective  elastic  modulus  much  greater 
tiumthatofice. 

The  ice  stress  is  determined  by  monitoring  the  radial  deformation  of  the  cylinder  with  a  vibrat¬ 
ing  wire  (Hawkes  and  Bailey  1973).  Savin’s  (1961)  stress-deformation  relations  for  cylindrical 
elastic  inclusions  in  elastic  and  viscoelastic  materials  are  used  to  calculate  the  applied  stresses 
from  the  radial  deformation  of  the  pup.  Three  sensors,  oriented  4S°  to  one  another,  are  used 
in  the  field  to  measure  the  magnitude  and  direction  of  the  principal  stresses  in  the  ice  sheet. 

The  pup  response  was  evaluated  in  the  laboratory  by  freezing  the  sensor  into  a  block  of  ice 
and  apidying  a  known  uniaxial  load.  Measured  stresses  were  pnerally  within  10%  of  the  applied 
stress  for  loads  up  to  2.1  MPa  and  ice  block  strains  less  than  0.25%.  In  long-term  creep  tests  where 
block  strains  of  4%  were  obtained,  measured  stresses  were  within  20%  of  the  applied  stress.  The 
pup  response  was  not  significantly  affected  by  variations  in  the  ice  modulus,  creep  and  differ¬ 
ential  thermal  expansion  between  the  ice  and  pup.  The  sensor  also  had  a  low  temperature  sen¬ 
sitivity. 

Hawkes  (1969b)  used  a  rigid  photoelastic  stressmeter  to  measure  stresses  in  blocks  of  frozen 
sand.  The  stressmeter  consisted  of  a  3.81-cm4ong  glass  cylinder  having  an  outside  diameter  of 
3.18  cm  and  a  wall  thickness  of  1 .27  cm.  hi  this  technique  stresses  are  determined  from  isochro- 
inatic  frinp  patterns  when  the  sensor  is  viewed  between  crossed  polarizers.  As  these  meters  had 
been  successfiilly  used  in  rock  and  concrete  under  elastic  loads,  the  objective  of  Hawkes’  experi¬ 
ment  was  to  evaluate  the  reqronse  of  the  sensor  in  materials  undergoing  creep  deformation.  In 
his  tests  the  sensor  was  inserted  into  blocks  of  frozen  sand  under  constant  load.  Measured  stresses 
were  found  to  be  within  10%  of  the  applied  stress  up  to  block  strains  of  5%.  Stiff  cylindrical 
sensors  have  also  been  tested  in  nonlinear  viscoek  -'ic  materials  under  both  uniaxial  and  biaxial 
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loads  (Hawkes  l%9b,  Skilton  1971 ,  Buswell  et  al.  1975).  Measured  stresses  were  within  2  to  10% 
of  the  applied  stress  for  both  short-term  and  long-term  loading. 

In  summary,  previous  investigations  have  shown  that  there  are  two  suitable  stress  sensor  designs 
for  measuring  stresses  in  ice:  a  stiff  cylindrical  sensor  having  an  effective  modulus  much  greater 
than  ice,  and  a  thin,  wide  sensor,  preferably  having  an  effective  modulus  close  to  that  of  ice. 

BIAXIAL  ICE  STRESS  SENSOR 

Hie  remainder  of  this  report  deals  with  the  testing  and  evaluation  of  a  stiff  cylindrical  sensor, 
described  by  Johnson  and  Cox  (1982),  that  is  used  to  measure  ice  stresses  in  a  biaxial  stress  field. 

It  is  an  extension  of  the  work  on  uniaxial  cylindrical  sensors  conducted  by  Johnson  and  Cox 
(1980).  During  the  preparation  of  this  report  we  learned  that  in  1975,  Ivor  Hawkes  of  IRAD  Gage 
suggested  to  Sun  Oil  that  a  stiff  cylindrical  sensor,  equipped  with  three  vibrating  wires,  could  be 
used  to  measure  the  principal  stresses  in  the  plane  of  an  ice  sheet  (Hawkes  1975),  As  Sun  Oil  did 
not  pursue  Hawkes’  suggestion,  the  development  and  testing  of  such  a  sensor  was  not  carried  out 
until  this  study. 

The  sensor  considered  in  this  investigation  consists  of  a  stiff  cylinder  made  of  steel  (Fig.  1  and 
2).  It  is  20.3  cm  long,  5.7  cm  in  diameter  and  it  has  a  wall  thickness  of  1 .6  cm.  The  ends  of  the 
sensor  are  threaded  such  that  a  rounded  end  cap  can  be  attached  to  the  lower  end  of  the  sensor. 
Extension  rods  can  also  be  screwed  to  the  top  of  the  sensor  to  position  the  sensing  portion  of  the 
gauge  at  any  desired  depth  in  the  ice  sheet. 

Principal  ice  stresses  normal  to  the  axis  of  the  gauge  are  determined  by  measuring  the  radial 
deformation  of  the  cylinder  wall  in  three  directions.  This  is  accomplished  by  use  of  vibrating  wire 
technology  advanced  by  IRAD  Gage  (Hawkes  and  Bailey  1973).  Three  tensioned  wires  are  set 
120^  from  each  other  across  the  cylinder  diameter  (Fig.  2).  The  diametral  deformation  of  the 
gauge  in  these  three  directions  is  determined  by  plucking  each  wire  with  a  magnet/coil  assembly 


Figure  1.  Biaxial  ice  stress  sensor. 
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Figure  2.  Schematic  of  biaxial  ice  stress  sensor. 


and  measuring  the  resonant  frequency  of  the  vibrating  wires.  A  thermistor  is  also  placed  inside 
the  cylinder  to  measure  the  gauge  temperature.  Both  ends  of  the  sensor  are  sealed  to  protect  the 
wires  and  electronics  from  moisture.  The  sensor  was  fabricated  by  IRAD  Gage  in  Lebanon,  N.H. 

This  design  offers  several  advantages.  The  sensor  is  rugged  and  leakproof  and  it  can  be  easily 
installed  in  the  ice  using  conventional  ice  augering  equipment.  As  the  sensor  output  is  in  terms 
of  frequency,  it  is  not  affected  by  leakage  to  ground,  poor  contacts  and  long  lead  lengths.  The 
sensor  is  also  inexpensive  ($1700  for  the  prototype,  including  labor  and  materials). 


BIAXIAL  STRESS  SENSOR  THEORY 

The  measurement  of  stresses  in  an  ice  sheet  with  an  imbedded  sensor  requires  precise  knowl¬ 
edge  of  the  strain  in  the  sensor,  the  stress-strain  relationship  of  the  sensor  material  and  the  sen¬ 
sor’s  inclusion  factor  under  different  loading  conditions.  Fortunately,  we  know  the  modulus  of 
the  biaxial  steel  sensor  and  we  can  precisely  determine  the  gauge  deformation  using  vibrating 
wire  technology.  Analytical  solutions  are  also  available  that  describe  the  behavior  of  a  cylindrical 
inclusion  in  a  plate  under  loading. 

Since  we  are  generally  interested  in  compressive  stresses  in  an  ice  sheet,  compressive  displace¬ 
ments  and  stresses  are  taken  to  be  positive  in  this  report  as  is  often  done  in  rock  mechanics. 
Principal  stresses  are  designated  by  p  and  q.  The  major  principal  stress,  p,  is  the  larger  compres¬ 
sive  stress,  sudi  that  p  >  q.  All  angles  are  measured  clockwise  from  the  p  direction. 

Geofe  defomation 

The  diametral  deformation  of  the  gauge  is  determined  by  measuring  the  resonant  frequency  of 
each  of  the  three  vibrating  wires.  The  fundamental  frequency  of  each  wire  is  proportional  to  the 
strain  in  the  wire  and  is  related  to  the  wire  strain  by  (Halliday  and  Resnick  1970) 
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/=  natural  frequency  of  the  wire  (s“*) 
=  wire  length  (5.08  x  10"*  m) 
e  =  wire  strain 
=  wire  modulus  (207  GPa) 

=  wire  density  (7.83  x  10*  kg/m*). 

Equation  1  may  also  be  expressed  as 


e  =  kf^ 


or  Ee  =  ~k  A(/*) 


(2) 


and/g  is  the  initial  wire  frequency.  Since  the  radial  deformation  of  the  cylinder,  K,,  at  e^/2  is 
equal  to 


we  have 


«w 


(3) 


Hence,  the  radial  deformation  of  the  cylinder  can  be  expressed  in  terms  of  the  change  in  frequency 
of  the  vibrating  wires. 

For  our  gauge  and  vibrating  wire  iheter.we  can  measure  radial  displacements  as  small  as  5.0  x 
10~*  pm  (2.0  X  10"*  in.).  This  corresponds  to  a  sensor  resolution  of  about  20  kPa  (3  Ibf/in.*) 
when  it  is  embedded  in  ice.  For  the  improved  vibrating  wire  meters  now  under  development,  the 
gauge  resolution  will  be  increased  tenfold.  Despite  the  stiffness  of  the  puge,  it  is  still  very  sensi¬ 
tive  to  loading. 


Stteaee  aaodated  wiA  cylindiical  sensors 

The  stress-deformation  relationship  for  cylindrical  elastic  inclusions  in  elastic  and  viscoelastic 
materials  has  been  examined  both  analytically  and  experimentally.  Savin  (1961),  Berry  and  Fair- 
hurst  (1966),  Williams  (1973)  and  others  have  developed  analytical  solutions  for  elastic  materials. 
Experimental  tests  have  verified  that  the  analytical  solutions  accurately  describe  the  stress  distri¬ 
bution  in  an  elastic  plate  (Suzuki  1969,  Wilson  1961).  The  analytical  solutions  also  describe  the 
defamation  of  cylindrical  elastic  inclusions  in  viscoelastic  and  other  time-dependent  materials 
in  uniaxial  and  biaxial  loading  experiments  (Hawkes  1969a,  b,  Skiiton  1971 ,  Williams  1973,  Bus- 
weli  et  al.  1975,  Johnson  and  Cox  1980). 

The  stress  and  displacement  equations  used  in  this  investiption  to  describe  the  behavior  of 
the  biaxial  stress  sensor  and  surrounding  ice  are  based  on  the  work  of  Savin  (1961).  Savin  (1961) 
developed  a  set  of  analytical  equations  to  describe  the  behavior  of  an  elastic  ring  welded  in  an 
elastic  plate.  Even  though  ice  has  time-dependent  properties,  the  analytical  results  of  Berry  and 
Fairhurst  and  Ae  experimental  work  of  Hawkes,  Skiiton,  and  Buswell  indicate  Aat  Savin’s 
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Figures.  Pkn  view  of  cylindrical  sensor 
frozen  into  ice. 


equations  can  be  used  in  our  application.  Because  the  gauge  is  significantly  stiffer  than  the  ice, 
its  response  diould  not  be  affected  by  variations  or  changes  in  the  ice  modulus. 

Generally,  we  are  interested  in  measuring  in-plane  stresses  in  the  ice  sheet.  Consider  a  cylindri¬ 
cal  sensor  that  is  frozen  into  an  infinite  isotropic  ice  sheet  (Fig.  3).  The  sensor  is  oriented  normal 
to  the  plane  of  the  ice  sheet,  which  is  subjected  to  in-plane  principal  stresses  p  and  q.  The  sensor 
has  an  outer  radius i?,  and  an  iiuier  radius^,.  The  stress  (a,,  and  r^^jand  displacement  (K, 
and  Vg)  equations  for  the  sensor  (/f ,  <  r  <  i!])  in  polar  coordinates  are 


'  -  .  ^ 

^  f’  ^  *-"■  ■  ”  ^  SI 

♦  C,  (,r,  +  l)  ^  +C4  COS29  (7) 

^  (c'3(^.  +  3)^-C,  ^-C,(;r,.l)^+C„^)sin2fl.  (8) 

The  stress  and  dispbcement  equations  for  the  ice  sheet  (r  >  /fj)  are 


8 


^2  3 

1+C,  -j  sin2d 

(11) 

►■.=  *= 

*  ^  *. 
8#ij  2 

^  ^  t  C,(jr|*l)  cos2» 

(12) 

y  s  D 

8m, 

^2.. 

(13) 

The  coefficients  C,  throu^  C,  depend  on  the  sensor  geometry  and  the  material  properties  of  the 
sensor  and  ice  where 

(14) 

/|3  (1+^,) 

(15) 

■  <^■) 

m  ^2  1 

''3  ^  D 

,*(«■- 0(^-1) 

(16) 

(1+^,) 
^4  =  -2  Z) 

(17) 

C=  ^ 
n* 

(18) 

C,  =  2.2  ^2^  C, 
n* 

(19) 

(1+^,) 

C  =2  ' 

^  D 

[(|.l)(4.3,.)»»‘(l.2-.i)]»> 

(20) 

C,=2-2  -^—  (^-*)  (3n«-6«<+4n*-l)  +  n‘(n2-l) 

(21) 

C,  «  -2  +  2  -  l)  (4n‘  -7»i^  +4«2  -  1)  +  n^  (n*  -  1)  ^  +^, 


(22) 


bi  the  above  equations 


E 

2(1 +i») 

X  ■  ^  for  plane  stress 


JIT  «  (3  -  Av)  for  plane  strain 


E  is  Young’s  modulus  and  v  is  Poisson’s  ratio  where  the  subscripts  s  and  i  denote  the  material 
properties  of  the  sensor  and  ice  respectively,  d  is  the  angle  measured  clockwise  from  the  princi¬ 
pal  stress  direction  p. 

Before  we  use  these  equations  to  compute  the  ice  stress  from  the  radial  deformation  of  the 
gauge,  it  is  instructive  to  examine  the  general  behavior  of  the  gauge  and  ice  under  different  load¬ 
ing  cooditiotu.  Since  the  gauge  is  bonded  to  the  ice,  the  radial  stresses  in  the  gauge  and  ice  are 
equal  at  the  ioe*gauge  boundary.  Equations  4  or  9  can  therefore  be  used  to  evaluate  the  gauge 
response  to  variatitMis  in  the  ice  modulus,  bt  Figure  4  the  normalized  radial  stress  o,  at  the  ice- 
gauge  boundary  parallel  to  the  uniaxial  loading  direction  p  is  plotted  against  the  ice-gauge  modulus 
ratio,  E^jE^.  In  performing  diis  calculation,  we  assume  plane  stress  and 

d-O" 


and 

It  is  apparent  from  Figure  4  that  in  situations  udiere  the  gauge  is  significantly  stiffer  than  the  ice, 
the  gauge  leqwnse  is  not  appreciably  affected  by  variations  in  the  ice  modulus.  The  ratio  ajp 
mnains  constant  at  a  value  of  1 .5.  This  corre^Kxids  to  an  inclusion  factor  of  0.67.  For  the  biaxial 


HpareA.  Normalized  radial  stress  at 
ice-ptuge  boundary  parallel  to  uni¬ 
axial  loading  direction  versus  ice-gauge 
modulus  ratio. 
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ioe  ttre«  lentor  that  hat  a  modulus  of  about  200  GPa.  there  is  a  slight  variation  in  the  incluaon 
fiKtor  for  ke  moduli  ranging  from  0.7  to  10  GPa.  However,  this  can  be  partially  compensated 
for  by  selecting  an  ice  modulus  appropriate  for  tfw  type  of  ice  in  which  the  sensor  is  installed. 

It  should  also  be  noted  that  creep  of  the  ice  reduces  the  effective  ice-gauge  modulus  ratio  and 
tesuhs  in  only  small  changes  in  the  inclusion  factor. 

Frederking  (1980)  designed  a  much  softer  cylindrical  ice  stress  sensor  having  an  “effective” 
modulus  close  to  that  of  ice.  Even  though  this  design  minimizes  overloading  of  the  surrounding 
ice.  the  sensor  response  is  greatly  affected  by  variations  in  the  ice  modulus  and  creep  because  it 
is  relatively  soft.  Frederking  obtained  poor  re»ihs  in  sensor  verification  tests  where  the  ice  mod- 
lius  wu  not  known. 

in  addition  to  variations  in  the  ice  modulus,  we  are  also  concerned  about  overloading  the  sur¬ 
rounding  ke  when  we  design  ice  stress  sensors.  Figures  S  through  7  give  the  stress  field  in  the 
ke  surrounding  the  biaxial  ioe  stress  sensor  f(»  both  uniaxial  and  biaxial  loading  conditions.  In 
preparing  these  plots,  we  assumed  plane  stress  and  chose  the  ice  modulus  to  be  3.0  GPa. 


a.  Radial  itnu.  b.  Tangential  stress. 
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c.  Shear  stress. 


Figure  5.  Nonnelized  stress  distribution  In  die  ice  surrounding  the  biaxial  ice  stress  sensor  under 
uniaxial  loads. 
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Figure  6.  Normalized  radial,  tangential  and 
diear  stress  distribution  in  the  ice  surround¬ 
ing  the  biaxial  ice  stress  sensor  under  biax¬ 
ial  loads  fq/p  =  1). 


Figure  7.  Normalized  radial,  tangential  and  shear  stress  distribu¬ 
tion  in  the  ke  at  the  ke-gauge  boundary  under  uniaxial  and 
biaxial  loads. 


ExamiMtion  of  these  plots  reveals  the  postible  shortcomings  of  using  stiff  cylindrical  stress  sen¬ 
sors  in  ice.  Radial,  tangential  and  shear  stresses  in  the  ice  can  be  increased  by  50%  due  to  the 
presenoe  of  the  sensor.  This  an  lead  to  premature  failure  of  the  ice  around  the  gauge  and  errone- 
oa  stiMS  measiaements.  b  addition,  if  the  Ixnid  between  the  gauge  and  the  ice  breaks  during 
localized  fee  failure,  the  above  eqations  are  no  longer  appliable  since  the  ‘Velded”  boundary 
condition  would  be  violated.  Controlled  laboratory  tests  are  needed  to  establish  the  operational 
bmits  of  the  cyUndrical  stress  sensor. 

Dstannbation  of  fee  stresses 

Hm  magnitude  and  direction  of  the  principal  stresses  in  the  ice  are  determined  from  the 
measured  radial  deformation  of  the  sensor  in  three  directions,  b  the  biaxial  stress  sensor  the 
measurement  directions  are  120”  apart  and  from  eq  7  we  have  for  die  displacements  of  the  three 
wires 


■i4(p-rq)-fR(p.q)cos20, 


(24) 
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•  A (p  * q)  +  B(p - q)  cos202 


(25) 


V,^"A(p+q)  +  B(p-q)cos203 


Ri  Rl  Rc  Ri  R( 

[^3(Ar.o)^  ^ 


+60* 


63*61  +  120*.  (30) 

bi  the  above  equations,  6]  is  the  angle  measured  clockwise  from  principal  stress  direction  p  to  the 
measurement  direction  ;  2R^  is  the  average  length  of  the  vibrating  wires  determined  from  a 
gauge  calibration  test.  Ideally,  >4  and  B  should  be  calculated  for  each  wire;  however,  all  three 
wires  are  essentially  the  same  length.  All  remaining  variables  have  been  previously  defmed. 
Solving  for  p,  q  and  6,  we  obtain 


1  \Vr^-A{p*q) 

®‘’2~*‘‘|  B{p-q) 


because 


cos  (6)  *  cos  (-6). 

Equation  33  has  two  solutions.  If 

*  A(p  +  q)  +  fi(p-q)cos2(6|  +60*) 
then  6|  is  positive.  If 

-  y40»  +  q)  +  «(p-q)cos2(6,  +120*) 
then  6,  is  negative. 
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Figures.  Variation  of  A  with  ice  moduti. 


—Ice  — j 

_ \ _ L_ _ L _ I _ L 

lOr  lo'  10*  10*  10*  10* 

Ei(MPo) 

Figure  9.  Variation  of  B  with  ice  moduli. 


Hie  coefficients  i4  and  B  depend  upon  the  geometry  and  mechanical  properties  of  the  sensor 
and  the  mechanical  properties  of  the  surrounding  ice.  It  is  desirable  to  design  a  cylindrical  sensor 
such  that  A  and  B  are  not  greatly  affected  by  variations  in  the  ice  mechanical  properties.  For  the 
biaxial  ice  stress  sensor,  wu  found  to  be  equal  to  2.387  ±  0.01 5  cm.  A  and  B  for  the  biaxial 
ice  stress  sensor  are  plotted  against  Young’s  modulus  of  ice  for  different  Poisson’s  ratios  in  Figures 
8  and  9  reflectively.  In  the  elastic  modulus  range  for  ke,  A  and  B  diow  little  diange  with  E^. 
However,  A  and  B  can  change  significantly  widi  diange  in  v,  Poisson’s  ratio.  Since  relatively  little 
is  known  about  Poisson’s  ratio  for  ice  (Schwarz  and  Weeks  1977),  we  assumed  a  value  of  033  in 
this  investigation,  bi  the  event  the  ice  deforms  plastically  widi  v  *  O.S,  A  and  B  would  then  be 
in  error  by  about  15%.  Fw  uniaxial  and  biaxial  (qip  >>  1)  loading  conditions,  this  would  result  in 
measured  ice  stresses  that  are  about  15%  too  hi^. 

Gange  cdftnthMi 

The  biaxial  ice  stress  sensor  is  not  calibrated  in  ice.  It  is  calibrated  in  a  hydraulic  pressure  ceU 
(Rg.  10)  to  determine  the  initial  frequency  and  the  effective  length,  2R^,  of  each  vibrating  wire. 
The  gauge  is  radially  loaded  and  the  measured  deformation  of  the  gauge  is  compared  to  the  radial 
deformation  of  a  thidc-wall  cylinder.  The  radial  deformation  of  a  thick-wall  cylinder  under  ex¬ 
ternal  pressure  at  a  radius /t,  is  given  by 


E,iRl-Rl)  '•'J 

adiere  Rj  and  i?,  are  the  outer  and  irmer  radii  of  the  cylinder,  p  is  the  radial  hydraulic  pressure, 
and  and  a,  are  Young’s  modulus  and  Poisson’s  ratio  of  the  sensor  material.  R^  corresponds 


B(inVlbf ) 


mi 
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Figure  10.  Hydraulic  pressure  cell  used  to  calibrate  the  biaxial  ice  stress  sensor. 


to  the  effective  radius  of  the  vibrating  wire  and  from  eq  2  and  3  we  have 
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A(P). 


From  the  calibration  test  we  also  have  the  alibration  coefficient  k^,  which  is  determined  by  the 
least  squares  method  from  the  test  data 

(36 

'  A(f*) 

Combining  eq  34, 35  and  36  and  noting  that 

we  obtain 

After  determining  the  effective  wire  radius,/?^,  by  trial  and  error,  we  can  then  use  eq  35  to  accu 
rately  obtain  the  radial  deformation  of  the  gauge  given  the  change  in  vibration  frequency  of  each 
wire.  The  average  value  for  all  three  wires  is  also  used  in  eq  27  and  28  to  calculate  A  and  B. 


EVALUATION  OF  THE  BIAXIAL  ICE  STRESS  SENSOR 

We  conducted  controlled  laboratory  tests  to  evaluate  the  biaxial  ice  stress  sensor.  The  tests 
were  first  conducted  to  determine  the  temperature  sensitivity  of  the  gauge.  The  sensor  was  then 


froam  into  luge  ic«  blocks  and  placed  in  a  biaxial  loading  machine  to  study  the  response  of  the 
semor  under  different  loading  conditions.  We  also  examined  the  effects  of  differential  thermal 
expensioa  between  the  sensor  and  sunounding  ice,  and  long-term  sensor  drift. 

Teugenituw  aemitivity 

We  determined  the  temperature  sensitivity  of  the  gauge  by  placing  the  sensor  in  a  glycol  bath 
inside  an  environmental  chamber.  The  temperature  of  the  chamber  and  bath  were  varied  and  sen¬ 
sor  readings  were  taken  at  different  temperatures.  The  results  for  each  of  the  three  vibrating  wires 
in  the  gauge  ue  presented  in  Figure  1 1 . 


Figwe  11.  Variation  of  wire  period  with  temperature 
for  each  of  the  three  wires  in  the  biaxial  ice  stress  sensor. 

Bfexial  hMNiing  test  equipment 

We  evaluated  the  response  of  the  sensor  to  ice  stress  by  freezing  the  sensor  into  large  ice  blocks 
and  placing  the  blocks  in  a  hydraulic  biaxial  loading  machine. 

The  large  ice  blocks  used  in  these  tests  were  61  by  61  cm  square  and  18  cm  thick.  We  chose 
the  length  and  width  of  the  block  to  accommodate  the  entire  area  of  influence  of  the  sensor,  about 
ten  diameters  (Fig.  5a).  The  blocks  were  grown  in  a  coldroom  using  the  freezing  chamber  shown 
in  Figures  12  and  13.  The  freezing  chamber  consisted  of  a  lucite  box  and  coldplate  surrounded 
by  foam  insulation.  To  grow  an  ice  block,  the  box  was  first  filled  with  either  fresh  or  saline  water. 
An  aluminum  coldplate  was  then  placed  in  the  top  of  the  box  and  connected  to  a  glycol  constant- 
temperature  bath.  The  sides  of  the  box  and  coldplate  were  then  covered  with  insulation  and  glycol 
was  pumped  through  the  coldplate.  The  coldroom  was  maintained  at  about  O^C  to  minimize  any 
growth  or  melting  on  the  sides  of  the  ice  block.  Additional  details  on  using  this  technique  to  grow 
ice  can  be  found  in  Weeks  and  Cox  (1974). 

This  method  of  ice  growth  produces  columnar  ice  which  is  very  similar  to  that  found  in  natural 
ice  sheets.  The  grain  size  of  the  crystals  in  the  fresh  water  and  saline  ice  blocks  varied  between 
OJ  and  2.0  cm.  The  fresh  water  blocks  had  both  horizontal  and  vertical  c-axis  crystals,  while  the 
salifw  ice  blocks  had  predominantly  horizontal  c-axis  crystals.  The  c-axis  of  the  crystals  did  not 
show  any  preferred  alignment  in  the  horizontal  plane.  A  vertical  thin  section  from  a  fresh  water 
ice  block  is  shown  in  Figure  14.  The  saline  ice  blocks  had  an  average  salinity  of  about  S  ^/oo. 
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Figure  12.  Freezing  chamber  used  to  grow  ice  blocks  for  the  stress  sensor  verifica¬ 
tion  tests. 


1 


Figure  1 4.  Vertical  thin  section  from  fresh  water  ice  block 
(scale  pi  cm). 

Hie  sensor  was  frozen  into  an  ice  block  by  drilling  a  hole  through  the  block,  supporting  the 
sensor  in  the  hole,  and  filling  the  annulus  between  the  sensor  and  ice  with  fresh  ice  water.  The 
center  of  the  block  was  first  found  and  then  a  6.4-cm-dianieter  hole  was  drilled  through  the  block 
with  a  small  auger  and  electric  drill.  An  extension  rod  was  then  attadied  to  the  top  of  the  sensor 
and  the  sensor  was  suf^rted  in  the  hole  with  the  frame  shown  in  Figure  1 5.  A  rubber  gasket  was 
next  placed  around  the  bottom  of  the  sensor  to  close  off  the  aimulus  between  the  sensor  and  ice. 
After  the  gasket  was  in  place,  ice  water  was  poured  into  the  aimulus  to  freeze  the  sensor  in  place. 
The  ice  block  and  sensor  temperatures  were  about  >10'’C.  The  next  day  the  support  frame  and 
extension  rod  were  removed  and  the  ice  block  was  positioned  in  the  loading  machine.  After  24 
hours,  the  ice  water  had  frozen  and  the  sensor  and  ice  block  temperatures  were  again  at  -10°C. 
Thermal  strains  in  die  sensor  and  ice  had  also  relaxed  by  that  time. 

The  biaxial  loading  machine  used  in  the  stress  sensor  verification  tests  is  shown  in  Figures  1 6 
and  17.  The  madiine  consisted  of  two  0.4-MN'Cspacity  hydraulic  rams  supported  by  two  inde¬ 
pendent  I-beam  frames.  The  inside  ram  and  frame  rolled  on  casters  (Fig.  16)  to  minimize  shear 
stresses  on  the  block  during  biaxial  loading  as  well  as  to  compensate  for  any  lack  of  planar  square¬ 
ness  of  die  ice  blocks.  The  {riatens  consisted  of  aluminum  blocks  covered  with  sheets  of  Teflon 
(Rg.  1 7).  They  were  only  58  cm  wide  to  allow  for  about  3  cm  of  block  deformation  during  a  test. 
If  the  fdatens  were  the  same  width  as  the  ice  block  (61  cm),  the  corners  of  the  platens  would  come 
into  contact  during  biaxial  loading.  The  platens  were  also  free  to  rotate  in  the  vertical  plane  to 
corapeoMte  for  any  lack  of  end  squareness  of  the  ice  blocks. 

We  apfriied  loads  to  the  ice  blocks  using  the  rams  and  a  hydraulic  hand  pump.  We  used  control 
valves  to  direct  the  hydraulic  fluid  to  one  or  both  rams  and  a  hydraulic  dial  gauge  to  measure  the 


*• 


18 


Figure  15.  Support  frame  used  to  position  the  sensor  in  the  hole  while  it  was  frozen 
in  the  ice  block. 


Figure  16.  Biaxial  loading  machine  used  in  stress  sensor  verification  tests. 
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Figure  1 7.  Ice  block  and  sensor  in  biaxial  loading  machine. 


Figure  18.  Equipment  used  to  measure  block  strains  during  testing. 


load.  The  entire  loading  system  was  calibrated  with  a  load  cell  prior  to  any  testing.  Applied  uni¬ 
axial  stresses  on  the  ice  blocks  had  an  accuracy  of  about  20  kPa. 

In  a  few  tests,  block  strains  were  measured  with  Direct  Current  Distance  Transducers  (DCDTs). 
The  DCDTs  were  attached  to  the  platens  as  shown  in  Figure  18.  This  equipment  was  not  available 
for  all  of  the  tests. 

Biaxial  loading  test  results 

We  conducted  sensor  verification  tests  in  both  fresh  water  and  saline  ice  under  uniaxial  and  bi¬ 
axial  loading  conditions.  A  long-term  loading  test  was  also  performed  to  evaluate  the  effects  of 
ice  creep  on  the  sensor  response.  All  tests  were  made  at  -lO^C. 

Equations  31  through  33  were  used  to  calculate  the  magnitude  and  direction  of  the  principal 
stresses  in  the  ice  from  the  radial  deformation  of  the  gauge.  The  coefficients  >4  and  B  were  cal¬ 
culated  to  be  equal  to  1 .1 52x  10“*  ^  and  4.739x  10"*  *  m^/N  respectively.  They  were  determined 
from  eq  27  and  28,  assuming  an  ice  modulus  of  0.69  GPa,  an  ice  Poisson’s  ratio  of  0.33  and  an 
effective  sensor  wire  radius, of  2.387x  10"*  m.  The  effective  wire  radius  was  determined  from 
the  previously  described  gauge  calibration  test.  We  assumed  a  low  value  of  the  ice  modulus,  in  that 
we  expected  loading  rates  on  the  biaxial  loading  machine  to  be  relatively  small. 

Four  ice  blocks  were  used  in  the  loading  tests.  Blocks  1  and  2  were  fresh  water  ice  and  blocks 
3  and  4  were  saline  ice.  The  relative  position  of  the  sensor  to  the  loading  directions  in  each  of  the 
ice  blocks  is  shown  in  Figure  19.  In  Figure  19,  Oy^,  Og  and  are  the  applied  stresses  on  the  ice 
block  and  p  and  q  are  the  measured  principal  stresses.  We  chose  the  principal  stress,  p,  to  be  the 
maximum  compressive  principal  stress.  The  angle  6  is  measured  clockwise  from  p  to  the  direction 
of  wire  1  in  the  sensor.  In  blocks  3  and  4  we  chose  wire  I  in  the  sensor  to  be  parallel  to  the  o^  or 
A  loading  direction. 


Blocli  No.  3  BlOCfc  NO.  4 


Figure  19.  Position  of  sensor  relative  to  the  loading  directions  for 
each  of  the  four  ice  blocks  tested. 


The  loading  system  was  designed  to  minimize  die  development  of  shear  stresses  on  the  ice  block 
faces.  However,  the  ice  stress  readings  indicate  that  small  shear  stresses  (up  to  60  kPa)  were  present 
on  the  block  faces  during  the  biaxial  loading  tests.  As  the  applied  shear  stresses  could  not  be  inde¬ 
pendently  measured,  it  was  not  possible  to  compare  the  applied  shear  stresses  to  those  measured 
by  the  ice  stress  sensor.  In  determining  the  ice  stress  measurement  error,  we  therefore  assumed 
that  the  applied  shear  stresses  were  zero. 

The  first  ice  block  used  to  test  the  sensor  was  made  from  fresh  water.  The  sensor  was  frozen 
into  the  block  with  wire  1  in  the  sensor  oriented  30°  clockwise  from  the  A  loading  direction  (Fig. 
19).  In  this  first  test,  the  platens  were  fixed  and  not  free  to  rotate  to  compensate  for  the  lack  of 
end  squareness.  A  small  load  was  therefore  placed  on  the  block  in  the  A  direction  to  cause  the  ice 
to  creep  and  come  into  good  contact  with  the  platens.  After  about  0.24%  strain,  this  was  achieved 
and  the  block  was  loaded  in  the  A  direction  in  0.24  MPa  increments.  The  test  results  are  given  in 
Table  1  and  plotted  in  Figure  20.  At  the  conclusion  of  the  test  the  block  was  highly  fractured  and 
we  discarded  it. 

The  second  ice  block  was  also  made  of  fresh  water  ice.  The  sensor  was  frozen  into  the  block 
with  wire  1  in  the  sensor  oriented  5°  counterclockwise  from  the  A  loading  direction  (Fig.  19). 
Wheelbearing  grease  was  also  applied  to  the  platen  surfaces  to  further  reduce  any  shear  stresses  on 
the  sides  of  the  block.  We  used  this  block  to  examine  the  sensor  response  under  both  uniaxial  and 
biaxial  loads  (Oy^  »  Og)  as  well  as  to  changes  in  direction  of  the  applied  stress.  Since  the  platens 
were  still  not  free  to  rotate,  a  biaxial  load  was  first  applied  to  the  ice  to  cause  the  ice  to  creep  and 
come  into  good  contact  with  the  platens.  Once  reasonable  contact  was  obtained,  the  ice  was  loaded 
equally  in  the  ^4  and  B  directitms,  then  only  in  the  S  direction  and  finally  in  the  A  direction.  The 
results  are  given  in  Table  2  and  plotted  in  figure  21 .  Since  small  shear  stresses  can  be  present  on 
the  block  loading  faces  when  biaxially  loaded,  measured  loads  Oy^,  Og  and  ry^g  were  derived  from 
p,  q  and  B  using  Mcdu  circle  theory. 


Tabkl.  Compuiton  of  applied  and  measured  stiess  data  for  block  1 . 
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The  third  ice  block  was  made  of  saline  ice.  The  sensor  was  frozen  into  the  block  with  wire  1  in 
the  sensor  oriented  parallel  to  the  i4  loading  direction  (Fig.  19).  Wheelbearing  grease  was  not 
applied  to  the  platen  surfiices  as  we  regarded  it  as  being  too  messy.  However,  the  platens  were 
modified  and  now  were  free,  to  rotate  and  compensate  for  lack  of  end  squareness.  We  used  block 
3  to  examine  the  effects  of  ice  creep  on  the  sensor.  At  the  end  of  the  creep  test,  the  block  was 
uniaxially  loaded  until  it  could  no  longer  sustain  the  load.  The  test  results  are  given  in  Table  3 
and  plotted  in  Figure  22.  Due  to  ice  creep,  a  constant  load  could  not  be  maintained  on  the  ice 
block;  therefore  the  load  was  adjusted  to  its  initial  value  each  time  a  stress  reading  was  taken. 
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b.  48  to  96  hours. 

P^ure22.  Measured  vo'sus  applied  stress  for  block  3. 


The  fourth  and  last  ice  block  was  also  made  of  saline  ice.  The  sensor  was  frozen  into  the  block 
with  wire  1  in  the  sensor  oriented  parallel  to  the  A  loading  direction.  We  used  block  4  to  evaluate 
the  response  of  the  sensor  under  biaxial  loads  udiere  the  applied  loads  o^  and  were  not  equal. 
Such  tests  are  difficult  in  that  lateral  strairu  cause  the  ice  to  bow  out  in  Ae  minor  loading  direc¬ 
tion,  degrading  the  squareness  of  the  block.  The  test  results  are  given  in  Table  4  and  are  plotted  in 
Rgure  23.  Again,  measured  stresses  in  the  loading  directions  are  compared  to  the  applied  stresses 
as  small  shear  stresses  develop  on  the  platens  during  biaxial  loading. 


Tabled.  Cosnpuisonofap^ied  and  measured  strca  data  for  block  4. 
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Figure  23.  Measured  versus  applied  stress  for  block  4. 


Differantid  thcnml  ejqwiifioii 

We  also  used  block  4  to  study  the  effects  of  differential  thermal  expansion  between  the  ice  and 
gauge.  To  simulate  conditions  in  the  field  where  the  tensor  would  be  covered  with  an  insulated  in* 
strument  box,  the  block  was  insulated  on  the  top  and  bottom  with  foam.  The  temperatures  of  the 
sensor  and  block  were  then  varied  from  •■20"C  to  O^C  and  back  to  -lO^C  to  examine  the  re^nse 
of  die  gauge.  After  the  stress  measurements  were  corrected  for  changes  in  temperature,  the  stress 
readinfi  were  generally  within  or  equal  to  die  resolution  of  the  sensor,  20  kPa. 

Lmig4em  drift 

During  the  course  of  the  evaluation  study,  sensor  readings  were  obtained  at  20**C  to  examine 
the  stability  of  die  unloaded  gauge.  The  period  of  vibration  of  eadi  of  the  three  wires  in  the  gauge 


F^ure  24.  Variation  of  wVe  period  with  time  for  earii  of  the  three 
wires  in  the  biaxial  ice  stress  sensor. 
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it  plotted  againtt  time  in  Figure  24.  Over  a  period  of  about  200  days  wires  1  and  2  showed  a  slight 
decrease  in  the  wire  period,  while  wire  3  showed  a  significant  change. 

Dftcwioa  of  teat  veanite 

The  biaxial  ice  stress  sensor  hat  a  low  temperature  sensitivity  and  the  sensor  response  is  not 
affected  by  differential  thermal  expansion  between  the  ice  and  gauge. 

The  tensor  output  appears  to  vary  linearly  with  temperature.  If  the  sensor  were  in  ice,  the  ob¬ 
served  dianges  in  output  would  correspond  to  stresses  of  about  S  kPa/°C.  Relative  to  the  resolu¬ 
tion  of  the  gauge  (20  kPa),  the  temperature  sensitivity  is  small.  In  many  short-term  applications 
of  the  gauge,  temperature  corrections  would  not  be  needed.  However,  in  situations  where  large 
vartethms  in  ke  temperature  did  occur,  the  sensor  output  would  have  to  be  corrected  to  obtain  the 
hif^iest  possible  accuracy. 

Metge  et  al.  (197S)  postulated  that  the  response  of  a  steel  ice  stress  sensor  would  be  greatly 
affected  by  differential  thermal  ejqnnsion  between  the  ice  and  gauge.  However,  the  results  of  this 
investiption  and  those  of  Johnson  and  Cox  (1980)  do  not  support  this  hypothesis.  If  ice  were  an 
elastic  material,  differential  thermal  expansion  would  be  a  problem.  Because  ice  creeps  under  low 
stress,  localizied  thermal  stresses  in  the  ice  around  the  sensor  rapidly  relax  and  are  unable  to  build 
up  to  any  significant  value. 

The  results  of  the  loading  tests  indicate  that  the  sensor  responds  immediately  to  applied  loads. 

In  general,  measured  stresses  are  within  15%  of  the  applied  stresses  for  both  uniaxial  and  biaxial 
loading.  The  sensor  response  does  not  appear  to  be  alTected  by  ice  creep  and  when  the  applied 
stresses  are  removed,  the  measured  stresses  fall  close  to  zero.  Reliable  stress  measurements  are  also 
obtained  well  beyond  yielding  or  teilure  of  the  ice.  In  addition,  the  sensor  can  usually  determine 
the  direction  of  the  ap|died  stresses  to  within  S'*. 

Fsrt  of  the  observed  error  can  be  attributed  to  the  resolution  of  the  puge  (20  kPa)  and  the 
loading  system  (20  kPi).  Combined,  they  account  for  about  30  kRs  of  the  observed  differences 
in  the  applied  and  measured  stress.  These  differences  are  significant  at  low  stress  levek.  Errors 
are  also  introduced  by  poor  seating  between  the  tides  of  the  ice  block  and  platens,  and  shear  stresses 
on  the  block  sides  during  biaxial  loading.  The  block  3  and  block  4  results  show  that  stress  measure¬ 
ments  improve  significantly  when  the  applied  stresses  are  held  constant  on  the  block. 

If  it  were  possible  to  solve  the  seating  and  shear  stress  problems  associated  with  the  loading 
machine,  measured  stresses  would  probably  be  within  10%  of  the  applied  stresses.  This  postulate 
is  supported  by  the  results  from  block  3  where  a  uniaxial  load  was  maintained  on  the  ice  block. 
After  about  3  hours,  the  block  appears  to  be  properly  seated  and  the  difference  between  the  applied 
and  measured  stress  is  less  than  or  equal  to  the  combined  error  associated  with  the  resolution  of 
the  pugs  and  loading  system.  Since  a  biaxial  field  can  be  described  as  the  superposition  of  two 
normal  uniaxial  fields,  the  same  results  should  be  observed  in  a  well-desiped  biaxial  loading  test 
where  problems  associated  with  poor  seating,  bulging  of  the  ice  block  and  shear  stresses  have  been 
eliminated. 

ftriodic  measurements  under  no  load  reveal  dut  the  sensor  exhibits  long-term  drift.  With  time 
the  strain  in  the  vibrating  wires  inaeases,  resulting  in  a  deaease  in  the  wire  period.  According  to 
IRAD  Gap,  the  pup  fabricator,  this  behavior  is  caused  by  outward  displacement  of  the  clamps 
holdhig  dte  wires  in  the  sensor.  Work  on  borehole  stressmeters  indicates  that  the  problem  can  be 
elfaninated  by  heat-treating  the  sensor  after  fabrication.  This  was  not  done  for  the  prototype  sensor. 


CONCLUSIONS 

RelMle  ice  stress  measurements  can  be  obtained  by  measuring  the  diametral  deformation  of 
a  stiff  steel  cylinder  embedded  in  the  ice.  By  measuring  the  deformation  of  the  cylinder  in  three 
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directioils,  we  can  determine  both  the  magnitude  and  direction  of  the  principal  stresses  in  the  ice. 
An^tical  solutions  describing  the  behavior  of  an  elastic  ring  welded  in  an  elastic  plate  adequately 
pre^t  the  sensor’s  inclusion  factor  (stress  concentration  factor)  in  ice,  despite  the  fact  that  ice  is 
a  tiine<dependent  material.  Since  the  sensor  is  considerably  stiffer  than  the  ice,  its  deformation  is 
not  signifkantly  affected  by  variations  in  the  ice  elastic  modulus  and  by  nonelastic  behavior.  It 
is  not  necessary  to  calibrate  the  sensor  in  ice. 

Our  controlled  laboratory  experiments  to  evaluate  the  biaxial  ice  stress  sensor  indicate  that  the 
tensor  has  a  low  temperature  sensitivity  (S  kPa/"C)  and  is  not  significantly  affected  by  differential 
thermal  expansion  between  the  ice  and  the  gai^.  Loading  tests  on  fresh  water  and  saline  ice  blocks 
containing  an  embedded  sensor  show  that  the  sensor  has  a  resolution  of  20  kPa  and  an  accuracy  of 
better  than  15%  under  a  variety  of  both  uniaxial  and  biaxial  load  conditions.  When  allowances  are 
made  for  poor  seating  of  die  ice  blocks  in  the  loading  machine  and  shear  stresses  on  the  platens, 
test  results  suggest  that  the  sensor  accuracy  may  be  better  than  10%  of  the  applied  stress.  Principal 
strms  directions  can  be  resolved  to  within  about  S". 

The  cyhndrical  sensm  does  not  gready  overload  die  ice  and  can  accurately  measure  ice  stresses 
well  be3rond  ice  yielding  or  failure.  The  maximum  stress  riser  produced  by  the  presence  of  the 
sensor  in  the  ke  is  about  1.5. 

The  sensw  is  also  rugged,  leak-proof,  and  can  be  easily  installed  in  an  ice  dieet  with  conventional 
apgering  equipment. 
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